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Bulk viscosity may serve to damp sound waves in a system of neutrinos coupled to very light 
scalar particles, in the era after normal neutrino decoupling but before recombination. We calculate 
the bulk viscosity parameter in a minimal scheme involving the coupling of the two systems. We 
add some remarks on the bulk viscosity of a system of fully ionized hydrogen plus photons. 

PACS numbers: 14.60.Lm, 98.80.-k, 95.35+d 
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INTRODUCTION 

It has long been known that a simple form of bulk vis- 
cosity is a potential source of dissipation in mixtures of 
relativistic and nonrelativistic fluids. For the cosmolog- 
ically interesting case of photons and matter in the era 
just before recombination, where shear viscosity does in- 
deed strongly damp waves for wave-lengths up to some 
critical length, Weinberg has given the expression for 
the bulk viscosity parameter, 
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Here r 7 is a mean free time for photon in the medium 
and the electron (and proton) number density n is held 
fixed in the derivative. Since the derivative in Q is of 
the order (1/3 — 10~ 10 ) in the early universe, in view of 
the large photon/electron ratio, the effect is absolutely 
negligible in this system. 

For what comes later in this paper, we comment on 
the parameter, t 7 . The development that led to Q was 
based in part on the solution by Thomas 0] of a relativis- 
tic transport equation for photons. However, Thomas' 
treatment included only photon absorption and emission. 
When Compton scattering is the interaction mechanism, 
it appears that t 7 should be the mean free time for Comp- 
ton scattering weighted by an efficiency factor for energy 
transfer, that is to say, multiplied by a number of order 
of the number of collisions required for temperature equi- 
libration of a system which initially has different temper- 
atures for the matter and for the radiation. When the 
scattering mechanism is Compton scattering and the sys- 
tem is near the recombination temperature this is a large 
number indeed. 

The motivation for the present paper is the possible 
existence of a scalar field, <j), with mass smaller than neu- 
trino masses and coupled to neutrinos. Over time, the- 
orists have found a number of reasons to entertain such 
a possibility Q- ||. The most stringent bounds on the 
coupling strength, at least until recently, have come from 
supernovae |l3j. Recently Beacom, Bell and Dodlc- 
son demonstrated the rather surprising result that, 
staying within these bounds, processes v + v — ► 4> + 4> 



can remove most of the neutrinos from the universe as 
the temperature declines into the region in which the 
neutrinos become non-relativistic, a region that roughly 
corresponds to the recombination region itself. However 
Hannestad 0] pointed out that the acoustic oscillations 
of the fluid of interacting neutrinos and 4>'s in this picture 
could alter the CMB angular power spectrum in ways 
that conflict with the data 1 . 

Beginning from this observation, Raffelt and Hannes- 
tad |17| have given coupling constant limitations that a 
true free-streaming requirement would put on the v — <j) 
coupling. But more detailed calculations of the effect of 
v scattering, annihilation or decay on the CMB and the 
matter power spectrum appear to show that the land- 
scape of coupling constant limitations is considerably 
more complicated that dictated by the "free-streaming" 
criterion, and, indeed, that in some coupling schemes, 
parameters that give a "tightly coupled" v — <p fluid at 
the recombination temperature are compatible with the 
data 0. 

The first purpose of the present paper is to point out 
a possible role for bulk viscosity in controlling the size 
of neutrino density fluctuations under some conditions. 
Clearly the region of application would be in between the 
tightly coupled case and the free streaming case. The 
origin of this dissipation is qualitatively similar to that 
which underlies JTJ. We have two gases, one completely 
relativistic, the other in the relativistic-non-relativistic 
transitional region. Thus if we apply an adiabatic com- 
pression, they change temperature by different amounts. 
There is a time lag in the equilibration, and consequent 
creation of entropy. Now, however, the analog of the fac- 
tor that produced the square of the electron- photon ratio 
in Q need not give an extremely small number, since the 
numbers of i/'s and ^>'s are comparable. 

In principle, a formalism for the evolution of inhomo- 



1 In ref. lUl from 1987 it had already been noted that an inter- 
acting sea of neutrinos and majorons could form structure in the 
universe; of course the details of the interesting application of 
these ideas are now different, in view of all of the data that has 
come in since then. 
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geneities that followed the detailed evolution of the po- 
sition and momentum distributions of each kind of par- 
ticle would have no need at all to mention viscosity ex- 
plicitly. Indeed, shear viscosity is subsumed within the 
standard evolution codes for the baryon, electron, pho- 
ton systems, entering at the quadrupole level where non- 
isotropic stresses are taken into account. Bulk viscosity 
is not included, nor need it be, in the system of matter 
plus photons; but interacting neutrinos can be another 
story. 

There are two important ways in which our v — <f> sys- 
tem requires computational mechanics somewhat differ- 
ent from those necessary to derive the Weinberg formula: 

1) In contrast to the case with matter and photons, 
where there is a conserved number of protons and elec- 
trons, the neutrino system will be taken to be symmetri- 
cal between particles and antiparticles, and the annihila- 
tion process into </)'s is an integral part of the evolution. 
This makes the determination of the analog of the factor 
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case. 

2) Looking at the parameters for the problem at hand 
it appears that the frequency dependence of the bulk 
viscosity will play a role in the application. The con- 
sequences of this dependence have been important in the 
role of bulk viscosity in damping the gravitational radi- 
ation instabilities in rapidly rotating neutron stars |l9| - 
|2^ . This dependence, being essentially a non-local ef- 
fect, does not rigorously fit into the framework that is 
best suited to the analysis of a relativistic fluid, viz, ad- 
dition of an additional term in the energy-momentum 
tensor in the comoving system, 



ST^ = -5 i < j ((uj)V ■ U , 



(2) 



where U is the space component of the four- velocity. We 
believe that it is legitimate, however, to crank the for- 
malism forward using J5J, and then, at the point when 
Fourier transforms are introduced in position space, to 
evaluate C( w ) at = v s k where v s is the speed of sound 
in the medium. 

In one respect, however, our v — <p system is simpler 
than the matter-photon system; the energy equilibration 
in a single collision is relatively efficient. This is because 
we have the direct transformation, v + v — > 6 + <b? 



2 In models that couple light scalars to neutrinos, the neutrinos 
might be taken either as Dirac neutrinos, with possible couplings 
that preserve lepton number, leading to u + v — » (f> + <f>, or as 
Majorana neutrinos that violate lepton number, in which case 
we have v + v —*</> + cf>. The amplitudes for the two cases are 
different (in angular dependence, for example) but we can take 
an estimate of total rate, as a function of coupling, G, to be the 
same in both cases. The distinction will not enter further in this 
paper, except here and there in terminology. 



The formalism that we shall present is completely ap- 
plicable to the case of photons and ordinary matter and 
this may possibly be of interest in other astronomical 
contexts. We shall apply our method to the ionized- 
hydrogen system, with a result related to 



CALCULATION 

We consider an adiabatic compression specified by a 
value of — 5V/V, where we do not allow energy transfer 
between the two species, leading to respective tempera- 
ture changes for the species of 5T^^. Defining the indices 
70 t „ for the respective species as, 
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we take 
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(3) 
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where it is to be understood that we consistently take real 
parts to get linear perturbations in the physical quanti- 
ties. The time rate of change of the temperatures is thus, 



(5) 



The rate of heat transfer per unit volume between the 
two baths is of the order of, 



Q = I>0 



Tn 



(G) 



where To is the mean temperature of both baths, the 
energy density of the <j) particle bath, and T an effective 
collision rate for 0's in the medium. By "effective" we 
mean the rate for collisions of a with the v cloud, sup- 
plemented by rate for (j> + (f> — * 2v, and multiplied by 
an energy transfer efficiency factor which is of the order 
of, but less than, unity, (and model dependent as well). 
Adding the effect of heat transfer to the right hand side 
of © gives, 
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Subtracting, defining ST(t) = T^(t) — T v (t), and solving, 
we obtain, 



ST(t) = Re 



iw(70 - 7 !/ )T 6 



w + r^([4]-i + [4]- 1 )T - 



(8) 



In the calculation of work, J PdV, below, the volume 
differential will be given by [d/dt(SV)]dt = bVu sin(wt). 
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We will therefore retain only the sin(wf) terms in the 
ST's, since the cos(uit) terms do not contribute to the 
dissipation. We denote these parts and the associated 
pressure changes by ST S , SP S . 

From we then find that the sin(wt) term in the os- 
cillation of the temperature difference is divided between 
the two baths as, 



ST$(t) = 5T'(t)- 



6T»(t) = -ST s (t)- 



X<t>h 



X<t>h 



»i 
-v i 



(9) 



Next we calculate the induced pressure changes of the 
two clouds, 



SP$(t) = 6T$(t) 



6P°(t) = 6T*(t) 



dPA, 



dT 
dPu 



dT 



a„6TZ(t). 



(10) 



Taking the </> mass to be zero, we list the three relevant 
parameters for the <j) cloud, 

^ = , = 3a , 70 = 4/3 , (11) 

where we use units h = c = fcs = 1. For the neutrino 
cloud in the domain in which the mass plays a role, nu- 
merical calculations of the corresponding parameters are 
described in the appendix. 

We write the net pressure oscillation amplitude as, 
5 Pi + 5P£, obtaining, 



SP S = 5 PI + 5P° = ~ub(smLjt . 



(12) 



where 

c 



^ 2 + r 2 ^([ c r]- 1 + [cn- 1 ) 2 T - 2 l 4 
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(13) 



The amount of mechanical energy converted to thermal 
energy in one period is given by 
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dt2(u> sm(ujt)dV(t) = 2irujb 2 ( , (14) 



giving the time averaged dissipation per unit volume, 



E/V = CwV/2, 



(15) 



identifying ( as the bulk viscosity. 

Using eq.(5) of ref.0 and staying in the temperature 
region that is (nearly) non-relativistic for the f's we cal- 
culate the reaction rate for d> + d>^i/ + i/m the mixture 



as T = G i Tg{m u lT) 1 where G is the v — \ coupling con- 
stant, and the dimensionless function g(m u /T) is given 
in the appendix. Introducing k ss 3w, we write £ in the 
form, 



c 



where 



H 



1 + k 2 /h% 
10- 4 ( 7l ,-l/3)(-^-l/3)T 2 



(16) 



(17) 



g(m v /T ) ([c^]"^^]- 1 ) P\ 
The "cutoff" wave number ka is given as, 

k = 10- 3 G 4 T ^(^) , (18) 

where 

F © =103[5( ^ /To)]1/23_1/2 



(19) 



F and H are dimensionless functions of to^/To, plotted 
in figs. 1 and 2 for the range of values in which we are 
most interested. They are both of order unity in the 
semirelativistic region in which there may be significant 
effects of bulk viscosity. The calculation of the compo- 
nents of the above functions, a„, 7„, and g is discussed 
in the appendix. 




FIG. 1: The function H(rn u /T) which appears in the bulk 
viscosity result, (IB) - 

The damping rate for a plane sound wave in a medium 
with bulk viscosity, is given by eq.(2.56) of ref.Q], 



damp 



k 2 C 



2(P + P) 



(20) 



We divide by the frequency for a sound wave of wave- 
number fc, take v s fa 1/V3, and shift to appropriate units, 
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and absorption (not scattering) provides the mechanism 
for energy transfer, and in which the photon absorp- 
tion rate is independent of energy, by beginning from the 
Boltzmann equation for the photon distribution function 



d_ 



f = r- 1 [(i- 



■;/Tr, 



)f-e 



(24) 



FIG. 2: The function F(m v /T) which appears in the result 
for the wave- number cutoff ko 1181 1. 



obtaining, 



rdam P /(fcv s ) ~ -28 



10" 



[lMpcJ-W G 



,2eV 



where 



fro 



[1 Mpc] 



30.3 



G 
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(21) 



(22) 



From 112 111 and (12211 we see that there can be strong damp- 
ing in several oscillations when the system is at the re- 
combination temperature, for a neutrino mass that is of 
the order of 1 eV and coupling G of the order 10~ 6 . 



MATTER PLUS PHOTON SYSTEMS. 

All of our development is applicable to the system of 
ionized hydrogen considered by Weinberg, and we com- 
pare our result for this case to (QJ in the low-frequency 
limit, £(lj = 0). This will provide both a check on our 
pedestrian formalism, and possible application to ionized 
hydrogen systems at much higher densities and temper- 
atures than prevail in the recombination region in cos- 
mology. In (|13J) . replacing the subscript (f> by 7 and 
v by "mat" , we use the adiabatic indices for the mat- 
ter, 7 ma t = 5/3, a mat = 2n e and the specific heat, 
Cy 1 = 3n e where n e is the density of the protons and 
the electrons. For the photons we use the numbers from 
(|11[) but with a and Cy doubled, since there are two 
polarization states. We obtain, 



C(o) = 



r(4aT 3 + 3n e ) 2 



(23) 



with a = 7r 2 /15, the black-body constant in our units. 
The Weinberg formula Q yields the same result if we 
choose T = 4t~ 1 . We can check the relation between 
these parameters, in a system in which photon emission 



We estimate the rate of energy transfer Q and from it the 
coefficient r in © by inserting / = (exp[w/T 7 — l) -1 , 
expanding to first order in (T 7 — T m ), multiplying by u> 
and integrating d 3 k. We obtain, 



r = 



t^W(T ) J 



duj LO 



[exp(w/T ) - 1] 



3.83T" 1 . 



where, 



W(T ) 



exp(w/T ) 



1 



(25) 



(26) 



We can directly apply ij^ to determine bulk viscos- 
ity in regions in which it is dominated by the "free-free" 
photoemission and photoabsorption process, rather than 
by Compton scattering. This process is completely ef- 
ficient in energy transfer, so that our previous caveats 
with respect to an efficiency factor do not apply. The 
results will be most interesting in a system in which we 



have neither 



<< n~ 



We take the 



photo-absorption rate, r a b s from the Kramers formula, 
modified by the Gaunt factor correction [2^]; 



TabsM 



16ne 6 n 2 p phi 



3m e uj 3 



-^ /(2T) M^), W 



where K$ is the modified Bessel function. 

Defining the parameter to be used in this case for T in 
© by Tff, we have, in place of (|2*5|l . 



W(T ) J 
= 114. e 6 m~ 



dujT ahs (uj) 
3/2 T- 7/ V. 



exp(u/T ) - 1 



(28) 



Taking T / m e as the efficiency factor for energy transfer 
in the Compton process, we find that the bulk viscosity 
is determined by Tff, rather than Compton scattering 
in all regions except those in which n e << T 3 . In the 
latter regions, however it is much smaller than the shear 
viscosity by virtue of the n e dependence of l|23|l . 



DISCUSSION 

The principal result of this paper is the expression for 
the bulk viscosity of the interacting gas of v's and (j>'s. We 
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took one neutrino flavor, and an initial state with van- 
ishing chemical potential, i. e. equal numbers of ^'s and 
i/'s (or equal numbers of right-handed and left-handed 
neutrinos in the Majorana variation). 

For the tightly coupled case with, say, G w 10 -5 , we 
would find very little dissipation in our models. For cou- 
plings of G m 10~ 6 , and a neutrino mass of order of an 
eV we find relatively strong damping of a plane sound 
wave in the medium. We note the following about this 
regime: 

1) In the simplest coupling schemes, at least, there 
would be a negligible number of cp's produced prior to 
neutrino-matter decoupling. But for G w 10~ 6 there 
would be production at a rate that serves to create an 
equilibrium distribution during the period in which the 
neutrinos are still completely relativistic and prior to 
photon decoupling. 

2) In this domain, the shear viscosity and bulk vis- 
cosity will make comparable contributions to the damp- 
ing of a plane wave. However, when the shift is made 
to a proper cosmological calculation, in which the shear 
viscosity is implicit in the way the scattering rate en- 
ters the quadrupole terms, and the bulk velocity damps 
monopoles, it is less clear to us what the relative contri- 
butions of bulk viscosity to the actual signal will be. 

As we move into the region G < 10 ~ 6 the dissipation 
becomes greater, and it appears that our effects would be 
an essential element of the anal ysis of the perturbations. 
We note that the authors of ref [Tj| found that their free- 
streaming criterion ruled out couplings down to a value 
G«10- 7 . 

Finally, as will have been obvious, we made no attempt 
to include the complications that would arise from three 
flavors of neutrinos with a mass splitting scheme to agree 
with the data, let alone the possible complications of <p 
couplings that carry anything other than the unit matrix 
in flavor space. When couplings that allow v decays are 
included, the basic physics of our description remains un- 
changed, but the numerology concerning the interesting 
domain of G"s becomes very different. 

I thank Nicole Bell for discussions. This work was 
supported in part by NSF grant PHY-04559 18. 



neutrino (and anti-neutrino, if we use Dirac z/s) number 
in our periodic changes we need to introduce an oscillat- 
ing chemical potential perturbation Sp v as well as an os- 
cillating value for 6V/V. For the Dirac neutrino case the 
Sp s are the same for neutrino and antineutrino (in con- 
trast to an equilibrium case, where they would be equal 
and opposite if there were zero net lepton number.) We 
now have an energy equation, 



and a number conservation equation, 



SV 



1 



~dT 



ST 



dn v 
dp v 



Sp. 



(29) 



(30) 



where p v is to be evaluated at zero after differentiation. 
Defining, 



A 



1 



Pv + Pu dT 



1 dn 
ndf 



rl dn 

n dp, 



1 



dp v 



p + p u dp v 



we have 



Sp v = AST , 



(31) 



(32) 



dpu . dp v 
dT dp v 



\Pv + Pu] 1 , 



O-v = ~~Qrp' ( 33 ) 



The quantities above are calculated from the expressions, 



1 



Pu 



Pv 



p 2 dp 



p dp 



1 



I + e (E-^)/T ' 



E 



p dp 



P 2 /E 



1 + e {E-n„)/T 



(34) 



(35) 



(36) 



APPENDIX 

We need the temperature response 7^ — 1 to an adi- 
abatic change in volume, and the pressure response a v 
to this temperature response, all of this with the interac- 
tions between v and <f> that transfer energy between the 
baths turned off. A complication is that in the true equi- 
librium state, in the presence of the annihilation mech- 
anism, the number of v J s and P's (the latter, if we use 
Dirac v's) are not conserved. The average state, to which 
we apply our periodic disturbance, thus retains the value 
p = for neutrino chemical potentials. But to conserve 



where E — p 2 + m v . The parameters 7„ — 1 , a v are 
functions of the single variable, T jm v . In figs. 3 and 4. 
we show plots of the two parameters. 

Note that 7„ approaches 4/3 at high temperatures, as 
expected. We note a borderline inconsistency in our pro- 
cedures, above, in that we will trust interactions that 
equilibrate the neutrino distribution thermally, such as 
v — v scattering (from <fr exchange) to do their job, while 
turning off (schematically) the change in neutrino num- 
ber and transfer of energy between the two seas, although 
all of these effects are of the same order in G. 

For the functio n qf rn^/T) we take the expression given 
in eq(5) of ref. [lj], for the v annihilation rate, and 
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FIG. 3: The index 7„(m„/T). 
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